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Abstract 

We study, via bosonization, the Landau fixed point for the problem of in- 
teracting spinless fermions near the Fermi surface in dimensions higher than 
one. We rederive the bosonic representation of the Fermi operator and use it 
to find the general form of the fermion propagator for the Landau fixed point. 
Using a generalized Bogoliubov transformation we diagonalize exactly the 
bosonized hamiltonian for the fixed point and calculate the Fermion propaga- 
tor (and the quasiparticle residue) for isotropic interactions (independently of 
their strength). We reexamine two well known problems in this context: the 
screening of long range potentials and the Landau damping of gauge fields. 
We also discuss the origin of the Luttinger fixed point in one dimension in 
contrast with the Landau fixed point in higher dimensions. 
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I. INTRODUCTION 



During the last fifty years the Landau theory has been a paradigm used to explain 
the experimental behavior and electronic properties of quantum Fermi liquids jl],0]. Ini- 
tially the theory appeared as a phenomenological framework with a few parameters fixed 
by experiments. The presence of unknown parameters reflected, at that time, the lack of 
a microscopic theory. However, it was an extraordinary and necessary first step. Landau 
himself also established the route for the microscopic explanation for the validity of the 
theory. The Landau theory became the main tool for the study of the effects of correlations 
in electronic systems and its foundation was eventually established on microscopic grounds 
using field theoretic methods f||-[5|]. 

Although the main idea behind the Fermi liquid theory is quite simple, that is, the idea 
of a quasiparticle, its realization in terms of microscopic calculations is far from that. The 
idea is that when an electron interacts with other electrons it polarizes its vicinity and a 
cloud is formed around it (in a electronic system a hole is formed around the electron due 
to the electron-electron repulsion 0). In its motion the electron carries an extra inertia due 
to the existence of the other electrons. In more usual words, the quasiparticle is a dressed 
electron. In this picture the interaction between the electrons does not affect the intrinsic 
properties of the electron (that is, its quantum numbers) but only its dynamics. In field 
theoretic terms, the existence of a quasiparticle is related to the presence of an isolated 
singularity in the one-particle Green's function ||. This singularity produces a Dirac delta 
peak in the spectral function of the Green's function. In the non-interacting case all the 
weight of the spectral function (or the quasiparticle residue, Z F ) is in the peak [Z F = 1). 
The ground state in this case is a filled Fermi sea (due to Pauli's exclusion principle) with 
a sharp singularity at the Fermi momentum which defines the Fermi surface. When the 
interaction is turned on the strength of the peak is weakened (0 < Zp < 1) but it is still 
infinitely sharp. The rest of the spectral weight is carried by an incoherent background 
which plays no essential role in Fermi liquid theory. The most important consequence of the 
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presence of this sharp peak is that the ground state is still a filled Fermi sea with a weakened 
singularity. The existence of quasiparticles is thus directly related with the existence of a 
well defined Fermi surface. 

Fermi liquid theory has also been used as the starting point in order to understand the 
nature of condensed states of matter such as the superconducting state. The theory which 
explains the behavior of superconductivity of simple metals, the BCS theory, starts from 
the fact that the normal state is a quantum liquid described by the Landau theory || . The 
properties of the superconducting phase thus depend on the properties of the normal phase 
of the material. 

The observed unusual properties of the normal phase of the cuprates |7]|| posed the 
question of the existence of new states of condensed matter which are not described by 
the Landau theory of Fermi liquids. In particular, the phenomenology of the normal states 
appears to indicate the vanishing of the singularity at the Fermi surface (Zp = 0) ||. Many 
theoretical scenarios for such a non- Fermi liquid behavior have been proposed PJIOH- 

The apparent failure of the conventional Landau theory of the Fermi liquid in the context 
of the cuprates (and, perhaps, in more general strongly correlated systems) has motivated 
a wide search for alternative theoretical tools that, in principle, could handle the effects of 
strong correlations. In one space dimension, where the interactions are always strong as a 



result of the kinematical constraints, bosonization [[II] has emerged as the main theoretical 
tool. Recently the problem of bosonizing a dense Fermi system in arbitrary dimensions has 
been the focus of intense research. The main ideas were introduced originally by Luther 



12], rediscovered recently by Haldane [13] and developed in great detail by Houghton and 



Marston [14| and by us [[Uj. Here we will use the method of bosonization by coherent states 



that we have developed recently JU 

In this paper we investigate a class of fixed point Hamiltonians for interacting fermions 
which exhibit Landau-type behavior. The main motivation of this work is to see how does the 
conventional behavior predicted by Fermi liquid theory arises in the bosonization approach. 
It is important to reexamine this well known problem not just as a check on our methods but 



also since, unlike the conventional many-body perturbation theory approach to Fermi liquid 
theory, bosonization is not based on self consistent resummations of perturbation theory In 
principle, it should yield exact results for the low energy behavior of the system. Thus, this 
is a necessary step if these methods are to be applied to more interesting physical systems 
in which the Landau theory is believed to fail, such as the problem of a dense system of 
fermions coupled to dynamical gauge fields. Such systems are central for the understanding 
of some of the most novel approaches to the problem of High Temperature Superconductors 
16fl and to the compressible states of the Fractional Quantum Hall Effect |17|JT8| . 

In previous papers |15] we have examined the Landau fixed point. Here we show by an 



explicit calculation that systems of fermions (at finite density, relativistic or not, continuum 
or lattice) which interact via scalar potentials in the absence of nesting or gauge fields 
belong to this universality class. We investigate the properties of the operators that create 
the physical low-energy states of the system described by Fermi liquid theory. We construct 
a Hilbert space of states which represents the physical states close to the Fermi energy. 
The fixed point Hamiltonians contain only marginal operators acting on these states. We 
also characterize the relevant operators at these fixed points which are connected with low 
energy instabilities of the system. We will not consider here the problem of spin and magnetic 
excitations. 

We have shown in our earlier work that it is possible to bosonize an interacting spinless 
fermionic liquid, at long wavelengths, in terms of operators which create particle-hole pairs 



close to the Fermi surface ||15|| . Our results showed that the bosons of the theory are sound 
waves which propagate on the Fermi surface and have no resemblance to free non-relativistic 
bosons. Essentially these bosons are topologically constrained to the Fermi surface and 
therefore they propagate in a non-flat metric. The dynamics of these bosons is related to 
the elastic properties of the Fermi surface, that is, the Fermi surface sustains a surface 
tension when the quasiparticle residue is non-zero (Zp 7^ 0) ||15|| . When the surface tension 
vanishes [Zp = 0) a phase transition occurs at the Fermi surface and many properties of 
the system change abruptly. We further showed explicitly that the bosonized theory yields 



the correct thermodynamic properties of Fermi liquids. 

In this paper we begin by reviewing in some detail the bosonization procedure we have 



introduced before ||15|| . We use this approach to develop the generating functional for the 
bosonic fields in terms of coherent states which are coherent superposition of particle-hole 
pairs and represent the distortions of the Fermi surface. This representation for the generat- 
ing functional allows us to discuss the form of the fermion operator in terms of the bosons. 
Here the similarity with the bosonization in one dimensional systems becomes immediately 
clear. We also show that bosonization and non- Fermi liquid behavior (or rather, non-Landau 
behavior) are not one and the same thing. We argue that bosonization is a more general 
concept and that the non-Fermi liquid behavior of one- dimensional systems is a product of 
the smallness of available phase space which enhances the interactions. In one dimension 
the constraints imposed by the conservation laws couple the oscillations of the Fermi sur- 
face (in this case only two Fermi points) in a manner which will be discussed below. In 
dimensions higher the one the number of degrees of freedom (or Fermi points) is infinite and 
conservation laws alone are not enough to drive the system away from the Fermi liquid fixed 
point. 

We rederive an explicit formula for the fermion operator written in terms of bosons and 
we use it to obtain the correct free fermion propagator in the limit of long wavelengths. 
This calculation reveals many important aspects of the bosonization procedure we have 
been developing and confirms, once more, the usefulness of this method. 

Since the boson operator is a product of two particle operators in terms of the original 
fermions, we show that a two particle interaction (which is written as product of four par- 
ticle operators) can be written only in four different forms which are bilinears in bosonic 
operators. Assuming an isotropic interaction for the fermions (that is, no dependence of 
the interaction on the position of the Fermi surface) we diagonalize exactly the problem in 
the thermodynamic limit. We calculate explicitly the fermion propagator and the quasi- 
particle residue for any strength of the potential. We show that, for local interactions, in 
dimensions greater than one, Fermi liquid behavior is always expected. Recently, a pertur- 



bative approach in terms of the bosonic operators has been developed in order to show that 
bosonization is able to reproduce the expected behavior of the self-energy for Fermi liquids 
in two dimensions |19| . Moreover, our methods give us hints for the possibility of breakdown 
of the Fermi liquid theory due to singular interactions at the Fermi surface |f20| . 



We reexamined the problem of dynamical screening of the fermion-fermion interactions. 
Unlike the conventional RPA approach, we do not resum the bubble diagrams self consis- 
tently since bosonization is a non-perturbative approach. Hence, one has to work with the 
bare interaction and let the dynamics of the system decide. Thus, instead of assuming a 
"screening first" scenario, we approach the problem of dynamical screening by calculating 
directly the full fermion one-particle Green's function. We show that dynamical screening of 
long range interactions is due to the excitation of particle-hole fluctuations with momentum 
transfer tangent to the Fermi surface. We further show that assuming screening first, as 
it is usually done in the conventional approach to Fermi liquid theory leads to physically 
incorrect results for one dimensional systems. 

We also discuss the problem of screening of external probes in the bosonic language. 
We show that while external scalar potentials are always screened, external gauge fields are 
not screened but get Landau damping instead (this is the case of a non-superconducting 
material). The reason is well known, liquids do not screen transverse oscillations at low 
frequency. We also show how bosonization can explain the differences between the screening 
of scalar (longitudinal) and vector (gauge) fields and in particular we obtain expressions for 
the susceptibilities (response functions) in both cases. We rederive the well known fact that 
the RPA result |2j is exact in the limit of long wavelengths. 

The paper is organized as follows: in section II we review our bosonization procedure 
introducing some details which were not discussed in our previous works; using the coherent 
states defined by the boson annihilation operator we obtain the generating functional for the 
bosonic fields in section III; in section IV, we discuss the form of the fermionic operator in 
terms of the bosons and its relationship with the coherent state path integral developed in 
the previous section; with this machinery at hand, in section V we obtain the non-interacting 
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one-particle Green's function; in section VI we study what kind of interactions the bosonic 
hamiltonian can have if we start with a two-body interaction between the fermions, we show 
that we have only four kinds of terms which are possible in the bosonic language and we 
are able to classify them; in section VII using a generalized Bogoliubov transformation we 
diagonalize exactly the bosonized hamiltonian for fermions interacting via isotropic interac- 
tions in dimensions higher than one; in section VIII we calculate the one particle propagator 
and the quasiparticle residue as a function of the potential strength for Fermi liquids; in 
sections IX and X we obtain the well known results for the response functions for scalar 
and vector fields, respectively, and in section XI we discuss the differences between one and 
higher dimensions in the context of bosonization. Section XII contains our conclusions. 



II. BOSONIZATION 



The bosonization of a fermionic system is based on the algebra obeyed by the densities 



and currents This algebra is obtained in a restricted Hilbert space which contains the 
states close to the Fermi surface. The relevant operator, which generates all the states in 
this Hilbert space is defined as, 



2 2 

In general the commutation relations between the operators in (2.1) are written in terms 
of an expansion of operators. However, in the restricted Hilbert space, we expect that 
the substitution of the commutation relations by their expectation value on the state \FS) 
(representing the filled Fermi sea) will generate all the relevant dynamics of the interacting 
fermionic system. The commutation relation is obtained for the case of small momenta 
normal to the Fermi surface (there is no restriction for momenta tangent of the Fermi 
surface). We have shown that, in this Hilbert space, the commutation relation between the 



operators in (2.1) is written as [15| 



n r -, 



{k),n_- g/ (k>) = SgoS^q- ^(/i - eg). (2.2) 



We begin by denning a complete set of one particle states with spectrum eg which is used 
to build the full Hilbert space. The velocity of the particles is defined in the usual way, 

n = Vcy. (2.3) 
The Fermi surface is defined by the set of vectors {kp} which obey the relation, 

A* = % (2-4) 

where [i is the chemical potential of the system. 

Although the particle-hole operators (2.1) have almost bosonic character, the operators 
do not annihilate the reference state \FS). We need to normal order these operator relative 
to this state. Also, canonical bosonic commutation relations are only obeyed by suitably 
smeared operators at each Fermi point. We define creation and annihilation operators, 

a q ik F ) =J2$A(\k-k F \) (n q ik)Q{v^ F ■ +n^k)Q{-v %F ■ q}) (2.5) 

k 

and 

4(M = E $ a(I^- k F \) (n_f$e(tT fc ■ q) + nfim-v^ ■ q)) (2.6) 

k 

where Q(x) = 1(— 1) if x > 0(< 0) and $a(|/c — is a dimensionless smearing function 
which keeps the vectors k close to kp that is, 

Lim A ^ A (\k -k F \) = % f (2.7) 

where A can be viewed as a cut-off in momentum space. The idea is to construct spheres of 
radius A which cover all the states on the Fermi surface |]I3"|| . We parametrize each one of 
these spheres by the Fermi momentum kp, with q at the center of the spheres (alternatively, 
we could also have constructed pill boxes of height A and length D instead of spheres |L4| , 
the difference will be immaterial since all the physical quantities will not depend on the way 
we introduce the cut-off). This construction will give good results whenever the states in the 
problem have momentum close to kp with fluctuations of order q such that q « A << kp. 
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It is straightforward to see that, by construction, we have, 

a^k F ) I FS) = 0. (2.8) 

And the smeared operators a^(k F ) are found to obey the commutation relations, 

yfe), atXk'F)] = N A (k F )V | q.v %F \ S Mr faj + , (2.9) 

where N\(k F ) is a local density of states defined as follows (V is the volume of the system), 

Mk F ) = \ £ I $ a(|£ - k F \)\ 2 5(n - e s ). (2.10) 

Eqs. (2.8) and (2.9) show that these operators have bosonic character and generate the 
restricted Hilbert space of interest. 

The local density of states is a measure of the number of states per unit of energy per solid 
angle in the Fermi surface. In general, in the absence of Van Hove singularities, the local 
density of states is well behaved and independent of the cut-off. In this case we substitute 
it by its natural average which is the total density of states N(0), 

N(0) = ^J2S(^-e n ) (2.11) 

k 

divided by the solid angle on the Fermi surface = f dfl. Indeed, in our previous articles 
where we have studied the transport and thermodynamic properties of Fermi liquids we have 
used N A (k F ) = N(0)/S d P]. 

Although we could work directly with the operators defined above it is usual to rescale 
the operators in such way to absorb the density of states in its definition. We set, 

b^k F ) = {N K {k F )V\q.v %F \y 112 a Q ik F ) (2.12) 

which obey the usual bosonic algebra (independent of the cut-off), 

Uk F ), bt,(k' F )} = h F ,9 F (<W + k-?) ■ ( 2 - 13 ) 

These relations between bosonic operators, which have particle-hole character, will be the 
basis of our work. 
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III. COHERENT STATES AND GENERATING FUNCTIONAL 



Although we have a closed algebra and a reference state to work out the physics, we 
do not have yet a clear physical interpretation for the bosonic operators. We know from 
eq. (2.1) that they are related to particle-hole excitations. However, in order to develop more 
physical intuition about the excitations created by these operators, we use the coherent states 
associated with them. 

The coherent states are defined via a unitary operator, 

/ 



WD = exp 



-E 



N A (k F )Vq ■ v- k 



^k F )n-^(k F ) 



(3.1) 



\ k F ,q J 'AV'^V y "k F 

— » — » 

where U = U* implies (p-^(k F ) = (jf^{k F ). The operator U is a functional of the fields 
<Pq{k F ) which are defined on the Fermi surface. 



We can also rewrite (3.1) as, 



W]) = exp 



- E 



1 



N A (k F )V\q- v- kp 



(0 g -(fc F )n_g-(A; F ) - ^{k F )n^k F )) 



y k F ,q,q-vj: F >0 

and using the definitions (2.5) and (2.12) we have, 

/ / \ 1/2 



(3.2) 



WD = exp 



- E 



1 



N A (k F )V\q-v- k 



k F )bUk F ) - <f>*Jk F )b q ik F ] 



J 



(3.3) 



y k F ,q,q-v£ F >o 

The coherent state is defined as the evolution of the reference state via the operator U, 



m) = um)\FS). 

It is easy to show that this state is an eigenstate of the destruction operator, 

1/2 



(3.4) 



i 



<MM I [</>]>• 



(3.5) 



N A (k F )V\q-Vj: F \ 

This last result has a clear physical meaning. It means that the coherent state represents 
a deformation of the Fermi surface at the point k F in the direction of q due to a coherent 
(collective) superposition of particle-hole pairs. Therefore, the bosonic field 4>^(k F ) is a 
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measure of this deformation at that point. These bosonic fields are topologically constrained 
excitations which propagate on the Fermi surface. 

As usual with coherent states, we can show that they are not orthogonal to each other 
and that they are overcomplete. The overcompleteness of these states physically means 
that the bosons which propagate on the Fermi surface are wavepackets [JT5]. This is fully 
consistent with the standard picture of the Landau theory. 

The propagation of the many-body system in time, from time to time t, can be obtained 
by a calculation of the S-matrix with initial and final states of the kind considered above. 
The S-matrix is defined in the following way, 

{M,t\[0\,O) = Me^M) (3.6) 

where H is the Hamiltonian of the system. 

Using the closure relation for these states it is easy to show that it can also be written 
in terms of path integrals, 



F 1 



M 

D 2 [<f>]e iSW (3.7) 

[01 



where D 2 U1 - FF d^k^d^jk^r) 
where U [0\ - ll fcj?)T ,^ >0 ■ k n a {% f )v\^v %f \ and 



iS[cf>} = £ AT 7f Z^TZ Z i \l (Hk F )cpt(k F , 0) + <p# F )<f>$ F , tj) + f drL{(j>{T)) 

N a I I v in ■ 11~ I L£ x ' ' JO 

(3.8) 



is the classical action for the motion of the bosonic fields and 

m = Y. \ (*&)^ - 0K&)^) - IBWr)] (3.9) 
is the lagrangian density. The rescaled Hamiltonian H is defined as 

H W r)] = £ N A (k F )VW ■ % ,| (3.10) 

(k F , [<p\\k F 



J k F \ 

kp 
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In the path integral the boundary conditions are defined by the S"-matrix (3.6), 

— * — * 

(p$(k F , 0) = (p^k F ) 

<f0 F ,t) = tf$ F ). (3.11) 

As in any field theory we can also write a generating functional, Z, which is the trace of 
the S'-matrix. In terms of path integrals it is written as, 

D 2 [4>]e^ dt L ^F,t)4> 3 {k F ,t)) (3.12) 

where the Lagrangian is same as in (3.9) and (3.10). Now we have periodic boundary 
conditions due to the trace. As usual the Euler-Lagrange equations for the classical action 
will generate the semiclassical dynamics for the problem. We have shown that in case of 
a Fermi liquid the semiclassical dynamics is represented by the Landau equation of sound 
waves. Therefore, the bosonic fields represent these waves which are distortions of the 
Fermi surface evolving in time. Moreover, the same functional integral in imaginary time 
reproduces the thermodynamics of the Fermi liquids |T3 . 



As we will see in the next section these results, besides providing a new physical insight 
into the physics of fermionic systems (and an interpretation for the existence of bosons 
in a fermionic theory), are also powerful tools which will allow us to calculate correlation 
functions of interest. To this end, however, we need a dictionary to translate from the 
language of fermions to bosons. The aim of the next section is to provide it. 

IV. THE FERMION OPERATOR 

The connection between bosons and fermions is inspired from the bosonization methods 
in one-dimensional systems. Luther \T2] was the first to explore the analogy between one 



dimension and higher dimensions. Luther's idea was to define at each point of the Fermi 
surface a one dimensional system using the radial directions. However, Luther worked only 
with non-interacting fermions. Here, we present an argument paralleling Luther's work. 
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The fermion operator ip{r) is written in momentum space as, 

Inspired by the results in one dimension systems we follow Luther's construction of the 
fermionic operator \I>(r, k F ) as follows, 

kp 

where, 

V(r,k F )=f(k F )e J ^\ (4.2) 

for some operators f(k F ) and J(f,k F ). The correct commutation relations for the 
ty(r, k F ) can be obtained by imposing the commutation relation [n$(k), f(k F )] = and 

— * — * 

[n$(k), J(f, kp)] = c — number. This choice seems to be the simplest possible. 

Since the operators n$(k) generate the Hilbert space of interest, the form of the fermion 
operator \l/(r, kp) will be defined by the commutation relation between the former and the 
latter. It is easy to show, using the fermionic algebra for the operators eg, that, 

[^(r),^n q ik)]=e-^m- (4-3) 
k 

From (4.2) we have, 

[*(?, k F ),n q ik'p)] = e-^*(r, k' F ) S^. (4.4) 

— * — * 

Using the choices for the commutation relations between the operators f(k F ) and J(r, k F ) 
with n$(k), we easily get, 

[*(f, k F ),n q ik F )] = [n q ik F ), J(r, k F )]V(r, k F ). (4.5) 

Therefore, the fermion operator will be well defined if we insure that, 

[n q ik F ),J(r,k F )] =e~ 1 ^ (4.6) 
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Due to the commutation relations (2.2) it is easy to see that the correct choice is, 

J(r, A*) = - E AT 7f Z^Z Z n ~^)- ( 4 - 7 ) 
q N A {k F )Vq-Vj: F 

Notice the similarity between the form (3.1) and (4.7). This is not a mere coincidence. 
The fermion operator is represented as a coherent state of bosons. Since, as we will see, the 
bosons diagonalize the problem of an interacting electronic system, the fermion operator is 
a non-perturbative object in the language of the bosons. Moreover, since the operator U is 
a functional of the fields we observe that the fermion operator \I/(r, k F ) is given by U for the 
following choice of the fields 

Uk'F)=e-^% F ^ F , (4.8) 

which means that while the fermions interact among themselves the bosons are free excita- 
tions in the Fermi surface. 

From this observation we conclude that the fermion propagator, 

K{r-P,t-t') = £ (FS\*\r,k F ,t)*(?,k F ,lf)\FS) (4.9) 
kp,k' p 

can be written as, 

K(r-r',t-t')= r(k F )f(k' F )(k F ,r,t\k' F ,P,0), (4.10) 

kp,k' p 

where \k F ,f,t) is the coherent state (3.4) with the prescription (4.8). The functions f(k) 
are written in terms of the total number of fermions in the system [jl2 , 22| . Their role is to 
insure that the fermion operators anticommute with each other. 

The one-particle Green's function can be obtained directly from the propagator in the 
case when the actual ground state of the system is the filled Fermi sea: 

G(f - P, t-t') = K(f- P, t - t')Q(t - If) - K(f- P, t - t')*Q(t' - t). (4.11) 

As we will show in the next section this is the case of the non-interacting electronic system 
and Fermi liquids. 
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V. THE GREEN'S FUNCTION FOR THE FREE SYSTEM 



The hamiltonian for a non-interacting electronic system is just the kinetic term, 

k 

Since the operators (2.1) generate the Hilbert space of interest, the form of the hamilto- 
nian (5.1) in this space will depend on the commutation relations between the hamiltonian 



and these operators. It is easy to show that [15 , 



Ho, b^{k F )\ = -\q- V£ F \b ? (k F ), (5.2) 
therefore the non-interacting hamiltonian can be written as, 

#o = E E W-n F \b^k F )b ? {k F ). (5.3) 

k F ?,<?■%, >o 

This result is only valid in the restricted Hilbert space of states which lie close to the Fermi 
surface. 

Therefore, the free system, which is composed by only the continuum of particle-hole 
pairs, is described by a set of independent harmonic modes at each point of the Fermi 
surface. These modes oscillate with arbitrary phase difference between them. It is an 
incoherent oscillation of the Fermi surface which vanishes on average. It means that, on 
average, the shape of the Fermi surface is kept constant, as expected. 

Given the hamiltonian (5.3) we can use the methods of the last section in order to 
calculate the non- interacting Green's function. The rescaled hamiltonian defined in (3.10) 
is written as, 

H[4>{t)\ = \q- v %F \cj>*fi F )Uk F ). (5.4) 
The functional integral we obtain is the same as for the harmonic oscillator and it is 



easily done |21 . 
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We have seen that the propagator is obtained in this language by using the following 
choice (4.8), 



' 21% J 



(5.6) 



where a (~ A -1 ) is a cut-off in the normal direction to the Fermi surface which we intro- 
duce in order to regularize the radial integrals (exactly as in one dimension). Making the 
substitution above we find, 



(k F ,f,t\k' F ,0,0) = exp I- ^ 



Q,Q-vg F >0 



3N A (k F )V, 



a\q.v c I 



— — — — -e 




\<l- v k F \ 


\W9 F \ 



a\q.v- | 
fc p 

V p 



Z °k F ,k'p c \v Uf \ j{q.r-\q-v %F \t) 



■ (5- 

1 1 ' V k F I 

Notice that for different points of the Fermi surface the integral diverges logarithmically 
due to the factor l/\q.vg \ which diverges in the the limit of q — > 0. Thus, we can write, 



(k F , f, t\k' F , 0, 0) = k F ^ k' F . 



(5.8) 



If k F = k' F we find the following integral, 

/ 



a\q-v t 



(k F , f, t\k' F , 0, 0) = <% F fc -, p exp 



E 



v 



(1 



i(q-r-\q-vt \t)\ & 



V 



19 • v k f 



. (5.9) 



In order to evaluate this integral we observe that the component of q normal to the Fermi 
surface dominates the behavior of the integral. It is natural, therefore, to split q = q~N + qr 
where q~N = (with n^ p = which is the normal component to the Fermi 

surface and the tangential component q~r which is defined as qr-n^ = 0. This choice can 
be viewed as a Fresnel construction of differential geometry. We have built a local reference 
frame on the Fermi surface which follows its local geometry. 
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Observe that the integrals over the tangent component can be easily evaluated. They 
have the form, 

/_ dq T 4- 2 e^ T = (S(x T ))) (5.10) 

which gives a negligible contribution except for xt = 0. We conclude, therefore, that the 
part tangent to the Fermi surface does not contribute to the long distance behavior of the 
Green's function and therefore can be neglected in this limit. However, it does contribute 
to the counting of states at the Fermi surface. Indeed, from (2.10), we see that the local 
density of states can be written as, 

N A $ F ) = ± I ^|$ A (M| 2 (5.11) 
VJ \ v k F \ 

where dS is the area element on the Fermi surface. However, the area on the surface is 
written as an integral over the tangential component of q alone. Therefore, the tangential 
part of the integral contributes to the density of states and not to the long distance behavior 
of the correlation functions. Thus we make the following substitution, 



J2 N A (k F )V\vj: F \ J dq N . (5.12) 
The same type of substitution gives the correct thermodynamic properties for these systems, 



namely, a linear specific heat proportional to the density of states |T5 
Substituting (5.12) in (5.9) we finally get, 



(k F , f, t\k' F , 0, 0) = S irfir exp ( - jH dq N (l - e^^V^ ^— j . (5.13) 
Now we use the integral, 

/ dx(l-e lxa ) =ln(l-i~), (5.14) 

and we finally conclude, 

(k F ,r,t\k' F ,0,0} = - T 5 fj 9 l ■ . (5.15) 

which is the correct asymptotic form for the non-interacting system in any number of di- 
mensions 1 12]. It represents a fermion moving with velocity v% , as expected. 
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VI. THE HAMILTONIAN 



We now discuss the possible forms of the hamiltonian which can appear in the problem. 
We assume that the hamiltonian is composed of two terms, the kinetic energy of the fermions, 
eq. (5.1), and a two-particle interaction, 

p,p',q 

The scattering processes described by the operators in (6.1) can be visualized by the pairing 
of the annihilation and creation operators which are depicted in Fig.l-a and Fig.l-b. In 
Fig.l-a, we have small momentum transfer q (forward scattering) between the pairs. In 
Fig.l-b the pairs exchange momentum p — p' can be of order of 2k F (backward scattering). 

The form of the kinetic energy in terms of the bosonic operators was discussed in the 
previous section and it is given by Eq.(5.3). In terms of the operators which generate the 
Hilbert space of interest, it is easy to see that the interaction can be written as, 

U = ^ E U^{q) n_ q ip)n q ip>). (6.2) 

p,p',q 

We can now parametrize the form of the hamiltonian in terms of processes involving 
particle-hole pairs close to the Fermi surface. From (2.5) we get four kind of terms (these 
are the only types of terms (scattering processes) which are present in the restricted Hilbert 
space), 

a\{k F )a^k' F ) (6.3) 
for q-v %p > 0, q-vg p > 0; 

a^k F )ai{k' F ) (6.4) 
for q-v %p < 0, q-vg p < 0; 



a^k^a^k'p) 
for q-v %p < 0, q-vg p > 0; 
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(6.5) 



a\{k F )a\{k> F ) (6.6) 
for q- v^ f >0, q- Vg p < 0. 

For a fixed direction of q the processes (6.3), (6.4) and (6.5), (6.6) are shown in Fig.2 and 
Fig. 3, respectively. 

Observe that the momentum q breaks the rotation symmetry 0(3) of the Fermi surface 
(for the isotropic system, of course) and introduces an axis in the problem. We can now 
divide the Fermi sea into "south" and "north" hemispheres. While the interactions (6.3) 
and (6.4) occur within one of the hemispheres the interactions (6.5) and (6.6) link different 
hemispheres of the Fermi sea. 

In Fermi liquid theory only the interactions (6.3) and (6.4) are present 0. The reason 
for that, in the absence of nesting and singular interactions, is very easy to understand. The 
scattering of a particle-hole pair initially located at vector kp to a new position k'p requires 
an amount of energy \q - (v% F — "%' F ) I- Since in perturbation theory the denominators in the 
perturbative expansion are written as on — q- (v% F — "%' F )' a ^ ^ ow ener §y5 the m ost important 
contribution, for fixed momentum transfer comes from regions of momentum space which 
v% ~ vg . This is only true, of course, if the matrix elements between the states considered 
are not singular. Indeed, it is worth noticing that the corrections to the specific heat due 
to the interaction between particle-hole pairs is calculated in this limit [^ |15"| , p3| and they 
agree very well with the experimental results. The presence of this process will give rise to 
terms of the form a) a in the hamiltonian. That is, the same form as the non-interacting 
system (except for being non-diagonal in the index k) [jnj. Therefore, these interactions 
never mix the creation and annihilation operators for the bosons. As a consequence, the 
only effect that these interactions can produce is a different phase of oscillation for the fields 
(see section IX below), or, in other words, to renormalize the Fermi velocity. 

However, the terms (6.5) and (6.6) which link the "south" and "north" hemispheres are 
more interesting. They will appear in the Hamiltonian as a combination of a) a) and a a. 
The consequence of these terms is dramatic. Besides changing the phase of oscillation, they 
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also mix creation and annihilation operators. It means that they introduce a renormaliza- 
tion of the amplitude of the fields by a kind of Bogoliubov transformation. If this process 
survives in the thermodynamic limit it can be shown to give rise to anomalous dimensions 
in the correlations functions. We notice here that this is the same kind of interaction which 
produces anomalous dimensions in one-dimension. 

Indeed, suppose we have the interacting part of the hamiltonian written for an one 
dimensional lattice as, 

H = UY.PnPn+l (6.7) 

n 

where p n is the charge density at the site n. The interacting term (6.7) gives rise to two 



different terms in the bosonized hamiltonian [11 



U{ P i{q) P i{-q)+P2(q)p2(-q)) (6.8) 

2U Pl (q)p 2 (-q), (6.9) 

where pi{q) and p 2 {q) are the Fourier component of the "right" and "left" movers, respec- 
tively (they are analogous of the "south" and "north" movers). 

The interaction term (6.8) is associated with the processes described by (6.3) and (6.4) 
which are responsible for forward scattering. The process (6.9) is related to the terms (6.5) 
and (6.6). The same effect will occur in higher dimensions, for analogous reasons, if we have 
a nested Fermi surface. Observe that in the case of nesting the vectors v% at one side of 
Fermi surface are the same and therefore processes which translate the particle-hole pair 
within one side of the Fermi surface cost no energy (\q- \ — Vg,} | = 0). The only other 
available process in the system transfer particle-hole pairs across the Fermi sea. They are 
described by (6.5) or (6.6) (or (6.9) in the one dimensional case which is a special case of 
nesting). This process will appear in the case of a nested Fermi surface even in the absence 
of singular matrix elements. We note here that the vanishing of the quasiparticle residue due 
to nesting was already obtained in numerical calculations in the two-dimensional Hubbard 



model close to half filling [24 . 
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Notice also that the processes in different hemispheres of the Fermi surface require an 
amount of energy (|<f • (v^ F — "%' F ) I) which, for fixed momentum transfer, does not vanish 
because v% ~ ~^y F - Therefore, these processes will be marginally irrelevant at low energy 
and will not contribute for the physics of the states close of the Fermi surface. Indeed, by 
a perturbative renormalization group approach to Fermi liquids, Shankar has shown 
explicitly the marginality of these processes. In other words, in normal conditions these 
processes will scale away in the thermodynamic limit in the bosonic form of the hamiltonian. 
We will show this result explicitly in the next section with the exact diagonalization of the 
bosonic hamiltonian. For these processes to become part of the bosonic hamiltonian special 
physical interactions are needed. For instance, even in perturbation theory, these processes 
will become important if the matrix elements between the states is singular at low energies 
and long wavelengths. It is clear that in this case the perturbative expansion fails. More 
importantly, it is possible to show that coupling a Fermi liquid to dynamical gauge fields 
generically leads to interactions of this type and, consequently, to a breakdown of Fermi 
liquid theory. We will study this problem using our methods in a forthcoming paper 

VII. DIAGONALIZATION OF THE BOSONIC HAMILTONIAN 

As it was explained in the previous section the hamiltonian for interacting fermions close 
to the Fermi surface can be rewritten in bosonic form as, 

H = H + ]T yjq- v %p q- v,: p N A (k F )N A (k F ,)U(q) (btfyb^) + b\{-k F )b (i {-k F ') + 

k F ,k F /,q 

bfikFM-kp') + b\{-k F ,)b\{k F )) (7.1) 

where the sum is restricted to vectors such that <f-"% F > and <f-"% ; > 0. We also have used 
the definition (2.12) and (2.5) and changed the sums to be restricted to just one hemisphere 
of the Fermi surface. The free term, Hq is already defined in (5.3) and we assume that the 
Fermi surface is round and the interaction is isotropic (these assumptions are made in order 
to simplify the calculations, however none of them are really important in what follows). 
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In order to make notation simpler we discretize the hamiltonian in such a way that the 
Fermi sphere is divided into N patches |nj with area Sd-iA.^ 1 such that the total area of 
the Fermi sphere is written as Sdkf^ 1 = Sd-iA. d ~ 1 N. From this definition the local density 
of states at the Fermi points is simply N\(kp) = N(0)/N where the overall density of states 

S k^~^ 

is N(0) = ^y VF • 111 this case the dimensionless coupling constant of the theory is given by, 

g(q) = N A (k F )U(q) = = (7.2) 

which vanishes in the thermodynamic limit for dimensions greater than one. This is the 



route for the stability of Fermi liquid theory as already discussed by Haldane pT3 HT5| , p5 |. 

From now on we consider the Fermi liquid problem in two dimensions, partly because of 
its relevance of the cuprates and because of its general interest [p|,|26| . To each point j in the 
Fermi circle there is an angle 9j = jj-j. Since in the sum in (7.1) is restricted to q ■ v^ f > 
the values of j the allowed values of j are — ^ + l<j<^ — 1. In the direction of q we 
have j = 0. In this notation we can rewrite, 



Q' % = cos(0j). (7.3) 
Moreover, we define the following notation for the bosonic operators, 

b^kj) = bj 

bgi-kj) = a-j. (7.4) 

Using the previous definitions it is easy to see that the hamiltonian (7.1) can be rewritten 
in the following form, 

H = Y,qvFH(q) (7.5) 

q 

where 



Tt = S{ \ b\bi + a\aA + g ^ V^sJ \b\bj + a\<ij + afij + b\a}A , (7.6) 



1,3 



where s» = cos(0j). In the last expression we have dropped the index q since the original 
hamiltonian is already diagonal in this index. Observe that the hamiltonian (7.6) describe 
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a set of coupled harmonic oscillators. In order to diagonalize this hamiltonian we define a 
generalized Bogoliubov transformation which mixes different points at the Fermi surface. 
We introduce two real orthogonal matrices M.u and Afu and two new bosonic operators $ 
and on such that, 



h = E (Muft + Una}) 
i 



(7.7) 



The commutation relation between the operators enforces that, 



kj = J2 (MuMji - MuNji) 
i 

[k, aj] = = E (MuAfji - NuMji) , 



a h a\ 



(7.8) 
Si k and 



where the new bosonic operators obey the usual algebra, f3i,/3 k 
[PhPk] = [ai,a k ] = 0. 

It is assumed that these new operators diagonalize the problem completely, that is, the 
hamiltonian is written now as, 



(7.9) 



where Si is the eigenfrequency. 

To find the equation which defines the above matrices we look at the commutation 

relation between the hamiltonian and the operators. Using the hamiltonian (7.6) and the 
definition (7.7) with the hamiltonian (7.9) we find, 



(7.10) 



Substituting (7.7) and taking the commutation relations with the operators ai and fli we 
find the following equations, 



(Si - Si) Ma = gY, V^s] (Mji + Mji) 

3 

(Si + Si ) Nu = -gJ2 V 1 ^ ( m ji + -AGO 



(7.11) 
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These equations define the Bogoliubov transformation. 

Notice that in principle we have two different kinds of solutions: the first one is related 
to the particle-hole continuum and it is defined for Si < 1 and the second is related to the 
collective modes and is obtained for Si > 1. All the solutions of the particle- hole continuum 
can be written in the form Si = cos(#/ + 5i) where Si is a correction of order ~ and they 
correspond to an angular shift of the positions of the points at the Fermi surface. In the 
thermodynamic limit (which is the limit of interest) they fill densely the Fermi surface and 
renormalize to the bare frequencies 0. The collective modes cannot be written as a cosine 
of some angle. They detach from the particle-hole continuum and represent high energy 
excitations which are separated from it by an energy gap which depends on the interaction. 
However, in principle, they are also present in the solution of the bosonic hamiltonian. In 
the case of the Landau theory there is only one collective mode (the zero sound for neutral 
systems or the plasmon for charged systems) which is born at the direction of the momentum 
transfer q, or in our notation, at j — 0. In some sense the Fermi liquid problem resembles an 
impurity problem in which the states of the particle-hole continuum get only a phase shift 
and a bound state (the collective mode) appears in the spectrum. We therefore divide our 
solution in two follows. 



Observe that in this case the equation (7.11) can only be inverted if we disappear with 
the singularity which appears for i — I. this is done be writing (7.11) as, 



A. Particle-hole continuum Si = s; 




il — 



Sl + Si' 



(7.12) 



where, 




(7.13) 



j 
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The first term in the expression for M. a represents the solution for i — I where we introduced 
an unknown factor Z t which must be evaluated. The second term must be understood as 
the principal value of the function, that is, it vanishes for i = I. We must comment here 
that we also need to avoid the term i = —I since s\ = s_j. However, this only introduces 
simple modifications in the algebra and does not affect the content of the results. It is easy 
to show that with the introduction of these terms the final result will be given in terms of 
symmetric and antisymmetric combinations of the matrices that we obtain. We will come 
back to this point latter in the paper. 

By substituting the expression (7.12) in (7.13) we end up with an equation for Zi which 
can be written as, 




(7.14) 



B. Collective mode Si = So 

In the case of the collective mode there is no divergence in the equation (7.11) and the 
solution is just, 

ka _ aV^Co 

•M-io — — 

Mo = -^°°, (7.15) 

OO + Si 

where Co is defined as in (7.13). By substituting these equations there we find, 

i^EA (7-16) 

j ^0 S j 

which defines S - The solution of this equation can be seen graphically in fig.4. For a finite 
number of points the solutions for the particle-hole continuum have a finite angular shift 
while the collective mode detaches from the continuum. In the thermodynamic limit the 
particle-hole continuum renormalizes to the bare frequencies while the collective mode gets 
a finite renormalization in the frequency. 
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The sum in (7.16) can be rewritten in a well known form if we go back to our original 
notation and take the thermodynamic limit (N — > oo). Using the definition of the beginning 
of the section we find, 



9 V s 2 - s 2 W {q) J-,/2 2tt s 2 - cos 2 (#) 



3 3 



Jo Ztc s — cos(o) 



where, 



J^-) = N(0)ff^^ (7.18) 
\qv F J J S d uj-q-v F 

is the RPA polarization function [|I|]. Therefore, equation (7.16) is nothing but the equation 
for the collective modes, as expected. In two dimensions the form of the polarization function 



is given by [15 



II(*) = N(0) (-y==0(l*l " 1) + ^==j ©(I - \s\) ~ lj , (7.19) 

where we have introduced a small imaginary part in the denominator in (7.18). Notice that 
we can solve equation (7.16) for So using (7.19) and (7.17). Using the real part of (7.19) 
(the principal value) one finds So(q) = ^^n$)u( ) wn ^ cn * s a l wa Y s greater than 1. 

The normalization coefficients, Cj, in (7.12) are already undefined. In order to calculate 
these coefficients we go back to the equations (7.8). Using the results (7.12) the first of those 
equations gives, 

/ ZC 2 — Z C 2 

Sij = 5 iij Z 2 C 2 g 2 s i + g 2 ^/s~s~ (a< + sj) 



s 2 - s 2 



V + 2C ° S ° I (7->0) 

Zw ° °)(sf-s 2 ) (S 2 -s 2 )(S 2 -s 2 ) ] - 



while the second gives, 



Zj Cj - ZjC 2 = „ 2Cfsi ZCqSq 

2 _ r>2 A-^l f 9 o\ I 9 9\ f nl 9\ / n9 9\ ' \ ' / 



r (a? - * 2 ) (a? - *?) (S 2 - s 2 ) (S 2 - s?) 
By substitution of (7.21) in (7.20) we find finally that 
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We can further simplify considerably our results in the thermodynamic limit. The only 
sum which appears in the final form of the matrices Mu and Mu is the one related to the 
polarization function (7.19). It is easy to show that, 

E ^L = »L( 2 '^L__ 2 „) = _ N , (7.23, 
y sj-s] 2tt V l Jo Sl -cos9 J V ' 

since the principal value of the integral in (7.23) vanishes for \si\ < 1 (see equation (7.19)). 
If we substitute these results in (7.12) one gets, 

1 Ng ^/s^si 



Mu = S i7 i + — 



N 1 + Ng si-Si 



Mu = — — ■ (7.24) 

N 1 + Ng Sl + Si { } 

Which is the final form of the solution for the particle-hole continuum. Observe that Mu 
has a diagonal term of order zero in N' 1 . If we have take into account the fact that 
si = s_; it is easy to show that this term has the form M. u = 4= (5^ + sgn(i — 1)8^1). This 
term corresponds to a rotation of the identity matrix in the bosonic representation. Any 
calculation of physical quantities with these matrix would give the same result as for the non- 
interacting system. In terms of the Landau theory the rotation of the bosonic eigenfunctions 
is essentially the Landau's adiabatic principle which states that the states of the interacting 
system and the non-interacting system can be linked adiabaticlly |l]J^]. The next term in 
J\Au is order iV" 1 and represents the dressing of the particle by the interaction, that is, the 
formation of the polarization cloud. Although this is an effect of order iV -1 in the matrix A4u 
it may be appreciable, since the number of points of the Fermi surface is N. We will show 
that this is the case for the calculation of the quasiparticle residue. Furthermore, observe 
that the matrix Mu has contributions of order iV -1 only and therefore this matrix is null for 
the non-interacting system in the thermodynamic limit. Moreover, the matrix Mu carries 
no resonant terms in the denominator of (7.24). This can be understood easily from the 
discussion of the preceeding section. From definition (7.7) we observe that while Mu links 

27 



points on the same hemisphere of the Fermi surface the matrix Mu links points in opposite 
hemispheres. What we are showing here is that, since we choose a regular interaction, in the 
thermodynamic limit the contribution from opposite points of the Fermi surface vanishes 
completely. All the physics of the Fermi liquids is defined in just one hemisphere of the 
Fermi surface. We will show in the last section of this article that in one dimension this is 
not the case because the Fermi surface has a finite number of points. 

Another interesting feature of the solution (7.24) is its dependence on the interaction. 
Observe that in terms of our original variables the coefficient of the correction iV -1 in the 
matrices can be written as j^^u^j- Thus, if we have a long range interaction, a Coulomb 
interaction for instance, this term is always finite and we can replace it for a effective 
interaction which is local at long wavelengths. In the next section we show that this term 
has close relationship with the quasiparticle residue. Moreover, as we will show in section 
IX this term represents the screening of long range interactions which arises naturally in the 
bosonic formalism. 

VIII. THE FERMION PROPAGATOR FOR A FERMI LIQUID IN TWO 

DIMENSIONS 

The fermion propagator can be evaluated by using equation (3.3) and the Bogoliubov 
transformation (7.7). It is a tedious, but trivial algebra to show that the fermion propagator 
for two points, i and j, at the same hemisphere of the Fermi surface can be written as, 



(i,r,t\j,0) =exp<Yl 




1 



MuMu + MuMu MjiMji + MjiNji 



while for two points in opposite sides of the Fermi surface we find, 
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Therefore our task is to calculate the sums in these two expressions. If we use the result 
(7.7) we find, 

M u Mjie iVFqSlt = 5 id e iVFqSlt + — - — ^HEL (Js^qt _ e i Sj v Fq t\ + 
i 1 -|- N g Si Sj 

( a V , — ^ s^'- 

+ {rnr s ) <«, - 55 ( «, - .,>• (8 - 3) 

Notice that the last term in the above expression has a double singularity at si = Sj and 
s/ = Sj. We can extract this singularity using the Poincare's theorem, 

V- ^ = V—^— (— —) + n 2 5 ( Sl - Si ) 5 (s, - s,) . (8.4) 

The double delta function can be rewritten in terms of the angles at the Fermi surface due 
to the definitions in (7.6). Indeed, it is trivial to show that, 

S ,« - «) = i (cos W - cos m =1^=1^=, ( 8 .5) 
where we have used 61 = If we substitute (8.4) in (8.3) we finally find, 

where we have included only the terms of order zero in iV -1 , that is, the only terms that 
survive in the thermodynamic limit. By the same token we can show that, 

£ UuU^^ = O (1) . (8.7) 

As we have explained before this result shows that processes which involve particle-hole 
pairs at opposite sides of the Fermi surface disappear in the thermodynamic limit. 

If we substitute (8.6) and (8.7) in the expression for the Green's function (8.1) we end 
up with the following expression, 

(8.8) 
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where (i, f, t\j, 0)o is the propagator for the non-interacting system calculated in (5.15). 
Observe again the form of the interaction coefficient which appears in the screening form 
as we have pointed out in previous sections. That is, long range potentials will not change 
the quasiparticle residue because, at long wavelengths the potential seen by the fermions 
is purely local. In this paper, in order to simplify the calculation we just assume a local 
interactions from the beggining, U (q) = U, however we must keep in mind that screening is 
a natural consequence of bosonization in dimensions higher than one. 

Applying the same methods we have used in section V we make the replacement, 

dq T N($)v F 



q\ A 2 



(8.9) 



and we are left with the integral, 



A 

dq N q N e- 9 -£ (l - e ^(^*-«i-«0) = A 2 - 2 « A 2 , (8.10) 

(1 — zA yVpt — Tij • T)j 

where we have disregarded the last term in (8.10) since we are looking at the asymptotic 
state of the system. 

Substituting (8.10) and (8.9) in (8.8) we find, 

(i,f,t\j,Q) = Z F (i,r,t\j,0) (8.11) 

where, 

is the quasiparticle residue. Observe that this form of the quasiparticle residue has clear 
non-perturbative character. In the limit of weak coupling, that is, N(0)U — > 0, we find, 

(8.13) 



which is the expected form of the RPA solution for local interactions [pU27|]. In the strong 
coupling limit, N(0)U — ► oo, it is easy to show that, 

^°' 78 ( 1 + ^)' (8 - 14) 
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which, once more, has a non-perturbative character. Observe that the quasiparticle residue 
never vanishes and it actually saturates. That is, the Fermi liquid state survives for any 
strength of the interaction. Also observe that the propagator for opposite hemispheres of 
the Fermis surface, (8.2), vanishes in the thermodynamic limit. 

We want to emphasize that the result (8.14) is not completely universal. The numerical 
factor in front of equation (8.14) can change due to geometric factors. We have used Hal- 



dane's construction |13 of a sphere with radius A at each point of the Fermi surface. In 
this construction the Fermi surface is locally flat. If we use another different construction, 



such as in ref. ||14|| , we would get a different numerical factor. This factor is irrelevant and 
lead to a small renormalization of the quasiparticle residue. But this is not the important 
point here, the most striking result is the survivor of the Fermi liquid state for any strength 
of the interaction potential in the case of isotropic, local, interactions. 



IX. SCREENING OF A SCALAR POTENTIAL 

In the previous section we have shown that screening appears naturally in the exact 
solution of the Landau fixed point if we use bosonization. In this section we reconsider 
the problem of screening of external probes by the interacting Fermi system in terms of the 
bosonized theory. This is a well understood problem in the framework of Fermi liquid theory 
and the results of our bosonized theory naturally agree with the Fermi liquid theory results. 

In this section we examine the behavior of the system under a scalar external field whose 
external hamiltonian is written as, 

H ex t = ^V ext (q,t)p(^ (9.1) 

where p(q) is the Fourier component of the charge density of the system which is written in 
terms of the operators (2.1) as, 

= (9-2) 

k 

Therefore, we consider external hamiltonians of the form, 
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H ex t = ^J2 V e*t(q,t)n q ik) (9.3) 

q,k 

which, in terms of bosons, is written as, 

— » 1 /2 

H ext = £ (9.4) 

q,k F ,v% F -q>0 \ / 

It will be useful later to observe that the electronic density is given in terms of the bosons 

as, 

p(q) = _ £ (N A (k F )V\q- v %f \) 1 ' 2 (b\{-k F ) + bf{k F )) . (9.5) 

k F ,v %F q>0 

Here we consider the effect of the potential (9.4) in a Fermi liquid in order to see what kind 
of effect it can cause. The hamiltonian of the system in this case is given by H = H FL + H ext 
where H FL is given in (7.1). 

The equation of motion for the fields is given by the saddle point equation for the action 
in (3.8), 

% 



-q- * Sr M% F ) - U(q) Y (\q- %||? • vp F \N A (k F )N A (k' F )) 

k' F 

(Mk'F) + 4%-U F j) - ( lq - V ^ A{kp) \ Vext{qi t) (9. 6) 



which, after a Fourier transform, gives, 



q-v %F \N K {k F )v) 1/2 4>^k F ,u) = N A (k F )\V ext (q,u) 

u q ' v % F i 

+U(q) E (I? • vg F \N A (k' F )V) 1/2 (Uk'F,oo) + F ,uj) \ . (9.7) 

k' F ) 

Summing over k F on both sides of (9.7) and using (9.5) we find, 

(p(q,u))=n(q,u)(V ext (q,u) + U(q)(p(q,u))) (9.8) 

where, 



m^)=EN A (k F ) t /_/'^ (9.9) 

k F 



^-q- v k F 
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is the polarization function we have found in (7.18). 

The interpretation for (9.8) is very simple: the external potential produces a polarization 
of the fermionic gas which shields the interaction at long distances. Therefore, instead of 
the bare external potential V ext (q, uj) a new polarization potential appears, U \q) (p(q, u)) , 
and the effective potential which is felt by the fermions is given by, 

U eff (q,u) = V ext (q,u) + U(q)(p(q,u)). (9.10) 

The effective potential can be obtained from the bare potential using (9.8) and (9.10). 
Solving (9.8) for the density one gets, 

(p{q,w)) = - . 9.11 

1 - U{q)U(q,uj) 

Substituting (9.11) in (9.10) we finally find, 

U eff (q,u)= 9.12 

where e(q,u) is the dielectric function of the fermionic system which is given by, 

e(g» = l-tf(g)II(g». (9.13) 

The expressions (9.12) and (9.13) are the RPA results for the electronic gas Since we 

are dealing with the exact diagonalization of the fermionic system at long wavelengths, it is 
natural to recover the RPA approximation as the exact result in this limit since, as it is well 
known, its long wavelength limit saturates the sum rules. 

In the limit of interest, namely, low frequency and small momenta, the polarization 
function can be easily obtained. From (7.18) we have, 

na^m/a, . ~w (9 , 4) 

where we have included a small imaginary part for the frequency, Q is the solid angle and 
6 is the angle between the Fermi velocity and q. Observe that the polarization function is 
only a function of We can also rewrite (9.14) as, 

33 



m u) = mo) ( v I f^%r cosw ( jil - com)) 



(9.15) 



where V means principal value. 



In the limit of — << 1 we obtain 



U(q,u) = -N(0) ( 1 + ^—") + O ) , (9.16) 

V v F qJ \v 2 F q ) 



where 7^ is a numerical factor which depends on the spatial dimensionality, d, of the system 
(71 = 0,72 = 1/71,73 = 1/2). The second term in r.h.s. of (9.16) is the Landau damping 
term due to the decay of the bosonic modes into the particle-hole continuum. Observe that 
in one dimensional systems the bosonic modes never decay due fact that 71 = 0. It means 
that the true excitations of the one dimensional system are collective modes. In higher 
dimensions this is not true because there is always a particle hole continuum. In section XI 
we will see that this result has important consequences for the calculation of the Green's 
function. 

At this point it is tempting to follow the conventional approximation of Fermi liquid 
theory and to argue that the fermion-fermion interaction should also get screened just as 
much as an external probe is. However this is a delicate argument which is only justified 
by its success in explaining experiments. One should keep in mind that this is essentially 
a perturbative argument, motivated by a partial resummation of the perturbation theory 
series (namely the sum of all the bubble diagrams or RPA). In earlier sections of this paper, 
we showed that since the tangential fluctuations do not change the energy, they mainly give 
rise to a density of states. In the conventional "screening first" argument of Fermi liquid 
theory (in the RPA approximation ||) one finds an effective electron-electron interaction of 
the form 

■^"'' g-M-nit,) ' (9 ' 17) 

If we substitute the asymptotic form of the polarization II the RPA expression for the 
effective potential becomes, 
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This expression shows that, within RPA, screening is an effect caused by the density of states 
of the Fermi system. Notice, however, that this is exactly the same expression that appears 
in equation (8.8) when we diagonalized exactly the bosonic problem. Thus, bosonization in 
dimensions higher than one produces screening naturally and we do not have to put it by 
hand as in the usual approach in condensed matter physics. Furthermore, notice that our 
arguments in this section do not depend on the dimensionality of the system, therefore, we 
conclude that in any number of dimensions the screening of an external potential has RPA 
character. However, as we will show at the last section of this paper, in one dimension the 
fermion-fermion interaction is not screened. In particular, this result leads to the vanishing 
of the Green's function for the Coulomb potential in one dimension. Once more, if we have 
screened the potential by hand we would get a wrong result. 



X. COUPLING TO GAUGE FIELDS 

In the last two sections we showed how screening arises in bosonization of fermionic 
systems with long range scalar interactions. The reason for that is that the fermionic sys- 
tem resembles a liquid which can sustain longitudinal oscillations. However, if transverse 
oscillations are present the physics of the system changes completely. 

Let us consider a system of fermions in D > 1 space dimensions. In one space dimension, 
all gauge fields are purely longitudinal and, hence, the interactions they mediate are equiv- 
alent to Coulomb-like interactions of the form discussed in the previous section. Suppose, 
exactly as in the previous section, that we couple the fermionic system via minimal coupling 
with an external vector field. The hamiltonian can be written as 



H = H + H F _ G (10.1) 
Hq is the non-interacting hamiltonian and 
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H F - G = 9^-1 ■ A-Mk) + E^T-V^^ + Hg ( 10 - 2 ) 

q,k q,q',k 

is the fermion-gauge part of the hamiltonian where m* is the mass of the fermions. The 
first and the second terms on the r.h.s. of (10.2) give the coupling between the gauge fields 
and the fermionic system (g is the coupling constant of the theory). The last term on the 
r.h.s of (10.2) is a pure gauge field term usually describing the energy of the gauge field or 
possible external current terms which generate the field in the system. 

Since we are dealing with small momentum transfer only, we notice that the dominant 

— # — # 

contribution from the second term comes from q = —q' (the terms with q = or q' = 
are not present in order to preserve charge neutrality in the system). Using X^n (/c) = Nf 
where Nf is the number of fermions we can rewrite (10.2) in terms of the bosons operators 
approximately as, 

H F -g ~9 E (V- v %p \N A (k F )v) 1/2 v %F _ ql ■ (A_ ? b\(k F ) + A^(k F 



+^EAl 2 + Hg 

q 



;io.3) 



The generating functional is obtained exactly as before and the equations of motion for 
the bosonic fields are simply, 



dt 



-q- n F (t>q{k F ) - g (\g- v£ F \N A (k F )v) 



1/2 ^ 



v k F -i ■ A <i- 



(10.4) 



And for the gauge fields we find, 

d 1 



E 



Aq<t) 



9E 

kp 



k F - 



\q-vt\N K (k F ) 



1/2 



V 



g 2 N f 

m*V 



Aq(t) 



+ 



<t>q{k F ,t) 



(10.5) 



where D^ 1 (^,9) is the bare propagator for the gauge field (which in general is a tensor 
with spatial components i,j = 1,2, • • -d where d is the number of spatial dimensions) and 
Jext(r,t) is some external current. Eqs.(10.4) and (10.5) are coupled. They can be solved 
by a Fourier transform. Solving for (10.4) we find, 
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Substituting this result in (10.5) one finds, 



'•j 



h t 

kp,l 



q ■ v 



kp 











V? q 
F 2 _ 


i 




i 



M<*>] 



9 2 n f 
m* 



A g iuj) .+ J ext {q,u) . (10.7) 



where n/ = Nf/V is the density of fermions. Eq.(10.7) can be rewritten in a more appealing 
form as, 



J ext (q,u) 



10.8) 



where, 



D„ 1 (w, <?)1 . . - W ^(0)4^,,- + e ^a(^; 
j 1 >j \ 



q ■ v 



kp 



kp 



to — q ■ v 



kp 



;io.o) 



is the effective propagator for the gauge fields ||. We have used that m* = riF/{v F N(0)). 
In general it is very difficult to calculate the correction to the bare propagator. However, 
for an isotropic system, due to the symmetry, only the diagonal terms survive. For long 
wavelengths (q « k F ) and small frequencies (u> << v F k F ) one finds, 



IT 1 (u, q)}..= [Do 1 (u, q)].- v 2 F g 2 (N(0) + U(q, u)) S M 



(10.10) 



where U(q,co) is defined in (7.18). 

In the limit of interest, : 
renormalized propagator as 



In the limit of interest, namely ^ << 1, we can use the result (10.16) and rewrite the 



'•J 



Do" 1 (a;,, 



■2 U r 

. . + %u G di j 

hj v F q 



(10.11) 



where 



2 ng*nfy d 2 . , 2 

w P g = m ; = *ld9 N(Q)v F 



(10.12) 
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is the plasma frequency associated with the oscillations of the electronic system due to the 
coupling to gauge field. 

Observe that, contrary to the scalar case of section IX, the zero frequency gauge fields 
are not affected by fluctuations of the Fermi system, that is, 

[D (0, q)] itj — [D (0, q)]^ . (10.13) 

The cancellation of the density of states term in eq. (10.11) implies that the only effect of 
electron correlations at low frequencies and small wavevectors is a damping (not screening) 
of the transverse gauge fields. This is the well known phenomenon of Landau damping of 
transverse gauge fields in metals. A non cancelling density of states would imply a gap in 
the spectrum of fluctuations of the transverse gauge fields and the expulsion of static gauge 
fields, namely, a Meissner effect. This is what happens in a superconducting state. 

Let us assume, for the sake of the argument, that the bare propagator has the following 
form, 

[A) ML = ,. 9 * [A) Ml, (io.i4) 



hj CO 2 — VqQ 2 



hJ 



a tensor whose 



where v g is the velocity of propagation of the gauge fields. D (uj,q) 
form depends on the choice of gauge and it has an analytic dependence in uo and q. The 
renormalized propagator has the form (up to analytic structure tensors) 

u) v G q +iu pGvFq 

The Landau damping introduces a new physical scale in the problem. Observe that the 
interaction is screened with a frequency dependent screening length of order, 

/ 2 \ 1/3 

i-M-l^l ^ 1/3 , (10.16) 

V U P G J 

which defines a new scale in the problem. In the strong coupling limit, that is, g — > oo, the 
plasma frequency is much larger than the characteristic frequencies of the system {uo p g » 
uj). In this limit the characteristic momentum of the system will be cut-off by the plasma 
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frequency, ^== > q and the asymptotic form of the propagator is dominated by the 
Landau term, namely, 

D{uJ ^ = vW-l^^ (10 - 17) 

V GQ luJ pG VFq 

In real space and time this form of the propagator implies that the gauge fields behaves 
diffusively. The same type of propagator is found in the RPA approach 0,BM . 



XI. THE ONE DIMENSIONAL CASE 

In the previous sections we have shown that in terms of response functions bosonization 
gives the same result of the RPA approximation. This is due to the fact that RPA fulfill 
all the sum rules at long wavelengths and low energies, which is the exactly the limit where 
bosonization can be applied. Moreover, RPA is valid for high densities which in our language 
means that kp is large compared with the fluctuations in the system. This would lead us to 
conclude, erroneously, that RPA and bosonization are one and the same thing. Actually, the 
RPA results are expected from the bosonization point of view. We have already shown that 
bosonization gives rise to screening naturally, something that you have to put by hand in the 
RPA approach. Moreover, RPA is only valid in the weak coupling limit, N(0)U << 1 |§, 
and, as we have shown before, this is not the case of bosonization. It is indeed well known 
that in one dimension RPA works fine for correlations functions while it cannot explain the 
absence of isolated singularities in the Green's function. In this section we try to explain 
the reason for this behavior by comparison with our results in higher dimensions. 

In one dimension all the calculations simplify enormously. The matrices Aiu and Mu 
reduce to numbers. The condition (7.8) can be rewritten in terms of a variable ((q) such 
that, 

M{q) = cosh(C(g)) 

Af(q) = sinh(C(g)). (11.1) 
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The other condition in (7.8) is automatically fulfilled. In the one dimensional case there is 
no particle hole continuum since the Fermi surface reduces to two points, that is, s± — 1. 
The collective mode equation (7.16) defines the eigenvalue for the collective mode, 



which is easily solved as, 



S (q) = Jl+2g(q) = yJl + N(0)U(q) (11.3) 

where we used (7.2) with N = 2. The last expression can be put in a more standard form if 
we use the notation of section VII and rewrite the frequency of oscillation of the collective 
mode as, 



E q = v F qJl + N(0)U(q) (11.4) 



which is the well known result for one- dimensional systems | 2§|| . The variable ((q) is defined 
by the solution (7.15), that is, 



tanh ((? = -tttt = ~^—7 = ~ I n - 5 

M (l) S o + 1 Jl + 2g(q) + 1 



which can be rewritten in a more standard form, 



tanh(2C(g)) = r |^y (11.6) 



which is the expected result [29] 



Observe that M is finite, contrary to higher dimensions. This is a result of the finite 
number of Fermi points. We will now see that this has deep consequences for the one particle 
propagator. Indeed, using equation (11.1) we find, 

(i,x,t\j,0) = (*,Mi,0) exp IE i 1 ~ Ml^ iqVFSoiq)t )\ , (11-7) 
at the same Fermi point (right movers, for instance) and, 

(z,x,t|j,0) = (11.8) 
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for opposite Fermi points. This is the well known result for the one dimensional Luttinger 



model |29] 



The first consequence of the finite value of M is the presence of an anomalous dimension. 
Indeed, suppose the potential is local. Then the only effect of the interaction on the spectrum 



is a renormalization of the Fermi velocity from Vp to vp = vp^jl + N(0)U. In this case the 
integral (11.7) is easily done (see (5.14)) and the result is, 

a ( a 2 \ 7 , 

(i,x,tj,0) = , - 11.9 

x — vpt + la \x z — [vpt — %ay J 

where 

7= 4±*ffi._ L (11 ,o) 

Observe that the above Green's function has an anomalous dimensions given by 7 and a 
branch cut in the spectrum instead of an isolated singularity. As we discussed before this 
is a result of the process (a forward scattering!) that links opposite sides of the Fermi 
surface, a process which is suppressed by the presence of the particle-hole continuum in 
higher dimensions (in the absence of nesting, singular interactions or gauge fields). 

But this is not the only difference between one and higher dimensions. In the form of 
the fermion propagator (11.7) the interactions are not screened. Suppose, for instance, that 
we have a Coulomb interaction, that is, U(q) = e 2 /q 2 . This is the case of the Schwinger 
model which was studied via bosonization by Kogut and Susskind [[32] . In real space 



the Fourier transform of the Coulomb potential gives rise to a linear potential and therefore 
to confinement. If we calculate the spectrum from (11.4) we get a massive relativist bosonic 
theory (the system has a gap at q = 0). It is easy to verify that the propagator in (11.7) 
has an infrared divergence at t 7^ and it vanishes in this limit, that is, the fermions decay 
and disappears completely from the spectrum, only the collective mode is left. However, as 
we have shown in the section IX the one dimensional fermion gas screens external probes as 
expected, that is, in the RPA form. Of course, if we have screened the potential first, as in 
the RPA approach, we would never get this amazing result. 
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Therefore, the existence of the Luttinger fixed point in one dimension and the presence 
of anomalous dimensions is just a result of the lack of phase space (or due to a finite number 
of points in the Fermi surface). 



XII. CONCLUDING REMARKS 



In our previous papers [15| we studied the transport and thermodynamic properties of 



Fermi liquids. In this paper we have studied the Landau fixed point of a Fermi liquid and the 
associated one-particle propagator using the method of bosonization in arbitrary dimensions. 

We have shown that it is possible to bosonize a theory of interacting fermions in the 
restricted Hilbert space of states close to the Fermi surface. The bosonization is based on 
the algebra of the particle-hole operators in this Hilbert space and in the introduction of 
bosonic coherent states which generate deformation of the Fermi surface. From the coherent 
states it is possible to define a generating functional which is a path integral over the histories 
of the Fermi surface. The fields which propagate on the Fermi surface are sound waves which 
can be viewed as coherent superposition of particle-hole pairs. 

We have shown that from the construction of the fermion operator via coherent states 
we obtain the correct one-particle propagator which represents a fermion moving with the 
Fermi velocity. We also discuss the terms which appear in the interacting hamiltonian and 
and their relevance for the fixed point. It is obtained that processes which involve particle- 
holes on one hemisphere of the Fermi surface are associate with the Landau fixed point 
and the presence of anomalous dimensions, and the absence of singularities in the Green's 
function is due to operators which mix different sides of the Fermi surface. These operators 
appear in one dimension and in systems with nested Fermi surface due to lack of phase 
space for scattering or in the presence of singular interactions (such as in the case of gauge 
interactions). It is important to stress here that we are including forward scattering only 
and therefore our results are valid for gaussian fixed points. 

We show that for the simple case of isotropic interactions the bosonic hamiltonian can be 
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diagonalized by a generalized Bogoliubov transformation which mixes different points of the 
Fermi surface. We obtain two different types of solutions which represent the particle-hole 
continuum and collective mode. We obtain the fermion propagator in the thermodynamic 
limit for local interactions. We show that in dimensions higher than one the fermion prop- 
agator has isolated singularities and the only difference between the non-interacting system 
and the interacting one is the presence of the quasiparticle residue. We evaluate the quasi- 
particle residue for any strength of the interaction (and thus showing the non-perturbative 
character of the bosonization approach and its difference from the perturbative approaches 
based on ressumation of diagrams, such as RPA) and we show that the quasiparticle residue 
is always finite, that is, there is no possibility of breakdown of Fermi liquid theory for lo- 
cal interactions, exactly as expected. Furthermore, our results agree with the perturbative 
ones in the limit of weak coupling. Moreover, dynamical screening is a natural result of the 
bosonization method in dimensions higher than one. That is, we obtain that long range 
interactions are screened and we do not have to assume it as it is usual in perturbative 
theories in condensed matter physics. 

We also study the problem of the response of the fermionic gas to scalar external probes 
and obtain the RPA result for screening. This confirms that the bosonization is getting the 
correct physics at long wavelengths and low energies since the RPA fulfill the sum rules in 
this limit. And we stress once more that it does not mean that RPA and bosonization are 
the same thing because bosonization is a non-perturbative method which is valid for any 
strength of the interaction. We also show that when fermions are coupled to gauge fields 
the RPA result is also valid, the interactions are not screened but there is Landau damping 
(except in one dimension). We calculate the form of the effective propagator for the gauge 
fields and find the expected form for RPA. 

Finally we compare the one dimensional problem, related to a Luttinger fixed point, 
with the Landau fixed point. We show that due to the finite number of Fermi points in one 
dimension the mixing of points across the Fermi surface exists and leads to the appearance 
of anomalous dimensions in the Green's function. In one dimension there is only a collective 
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excitation while in higher dimensions there is a particle-hole continuum which shares spectral 
weight with the collective mode. We also show that the perturbative approach of screening 
the fermion-fermion interaction by hand would lead to wrong results in one dimension. While 
the one dimensional fermionic system screens external probes in the usual RPA form it does 
not screen the fermion-fermion interaction. The bosonization method leads to the vanishing 
of the Green's function in the case of Coulomb interactions in one dimension. In higher 
dimensions, due to screening, the Coulomb interaction leads to an effective local interaction 
which does no harm to the Fermi liquid behavior. 
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FIGURES 



FIG. 1. a)Forward scattering; b)Backward scattering. 

FIG. 2. Processes described by eqs. (6.3) and (6.4); a) Forward scattering; b)Backward scat- 
tering. The dashed line indicates the evolution of the state. 

FIG. 3. Processes described by eqs. (6.5) and (6.6); a) Forward scattering; b)Backward scat- 
tering. The dashed line indicates the evolution of the state. 

FIG. 4. Solution of the equation (7.16). The dashed lines represent the eigenfrequencies of the 
modes for the non-interacting system, the circles show the eingenfrequencies for the interacting 
system, a) particle-hole continuum; b) collective mode. 
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